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I.  INTRODUCTION 


Often,  in  the  design  and  production  of  complex  system  such  as  computer 
network,  good  system  reliability  is  a  critical  requirement.  Network  reliability  is 
rarely  treated  explicitly  and  quantitatively,  however,  largely  because  of  problems 
of  computational  complexity  and  because  few  efficient  computer  programs  exist 
for  computing  network  reliability.  Only  in  the  last  few  years  have  research  and 
development  in  this  area  provided  usable  results. 

We  wish,  as  efficiently  as  possible,  to  determine  the  network  reliability  for 
specific  class  of  networks.  The  purpose  of  this  thesis  is  to  develop  new  reduction 
and  decomposition  techniques  for  computing  network  reliability  and  to  show  how 
computational  complexity  can  be  decreased  by  using  these  techniques. 

The  network  model  which  is  used  in  this  thesis  may  be  thought  of  as  a 
communication  network  with  directed  edges,  allowing  only  one  way  communication 
from  tail  vertex  to  head  vertex.  Additionally,  we  impose  the  restriction  that  the 
network  have  no  directed  cycles.  Such  a  model  may  be  appropriate  in  hierarchical 
communication  networks  such  as  command  networks.  Also,  such  networks  may 
arise  as  subproblems  while  computing  the  reliability  of  general  directed  networks, 
Agrawal  [1962]. 

Networks  can  be  divided  into  various  classes,  e.g.,  directed,  undirected,  cyclic, 
acyclic,  etc.  A  large  amount  of  literature  exists  on  each  class.  (See  Agrawal 
[1962],  Ball  and  Provan  [1961],  Satyanarayana  [1982],  Hagstrom  [1960]  and 
Satyanarayana  and  Prabahkar  [1978]  for  directed  networks;  see  Ball  and  Provan 
[1961],  Johnson  [1962],  Satyanarayana  [1982],  Satyanarayana  and  Chang  [1981], 
Satyanaryana  and  Wood  [1982],  Valiant  [1978],  and  Wood  [1982]  for  undirected 


networks;  see  Agrawal  [1982]  and  Satyanarayana  and  Prabhakar  [1978]  for  cyclic 
versus  acyclic  directed  networks.)  The  specific  model  we  will  consider  is  the 
directed  acyclic  network,  denoted  D{V,E)  where  V  is  the  set  of  vertices  and  E  is 
the  set  of  edges.  Each  edge  is  represented  by  an  ordered  pair  of  vertices  e  ~{u,v). 
The  ordering  implies  that  the  flow  of  communication  can  move  along  a  from  u  to  v 
only. 

Each  edge  s<  in  the  network  functions  with  probability  p*  and  fails  with 
probability  g4  =  1— p4.  All  edge  failures  are  assumed  to  be  independent  events.  For 
simplicity  in  the  initial  discussion,  vertices  are  assumed  to  be  completely  reliable. 
However,  the  extension  to  unreliable  vertices  will  be  handled  briefly. 

Vertex  u  can  send  communication  to  vertex  v  in  D  if  there  is  a  directed  path 
of  working  edges  from  u  to  v.  One  of  the  vertices  s  is  designated  as  the  “source" 
or  the  "root".  With  respect  to  the  source,  there  are  three  different  measures  of 
reliability  which  are  usually  studied  in  directed  networks: 

(1)  The  source-to-terminal  (ST)  reliability,  R,t(D),  is  the  probability  that  s  can 

send  communication  to  a  specified  vertex  t  in  D. 

(2)  The  source-to-all-terminal  (SAT)  reliability,  R,v(D),  is  the  probability  that  s  can 

send  communication  to  all  vertices  of  D. 

(3)  The  source-to-K-terminal  (SKT)  reliability,  RajriD),  is  the  probability  that  s  can 

send  communication  to  every  vertex  in  a  specified  set  K  (these  vertices  are 
referred  to  as  "K-vertiees")  with,  by  convention,  s  e/f. 

Of  course,  (1)  and  (2)  are  special  cases  of  (3).  We  consider  the  general  case  (3)  in 
this  thesis. 

Generally,  the  of  problem  computing  network  reliability  is  NP-hard.  Certain 
reduction  and  decomposition  schemes  exist  which  reduce  the  size  of  the  problem 
while  preserving  network  reliability,  e.g.,  Agrawal[1982],  Shogan[l978]  and  Wood 
[1982],  and,  for  special  classes  of  networks,  it  is  thereby  possible  to  compute 
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reliability  in  polynomial  time.  In  general,  however,  the  solution  time  grows 
exponentially  in  the  size  of  the  network.  In  this  thesis,  we  develop  a  method  for 
computing  reliability  which,  in  general,  is  exponential  in  the  size  of  D.  but  which  is 
polynomial  when  the  number  of  non-K-vertices,  i.e.,  .  |  V-K | ,  is  fixed. 

The  rest  of  this  thesis  is  outlined  as  follows. 

Chapter  2  contains  a  necessary  definitions  and  graph  theoretic  concepts. 

Chapter  3  develops  reduction  and  decomposition  techniques  to  reduce  the  size 
of  a  network  in  order  to  decrease  computation  time.  Some  of  these  techniques  are 
new. 

Chapter  4  shows  how  to  compute  network  reliability  exactly  when  the  number 
of  non-K-vertices  is  moderate.  A  truncated  version  of  this  technique  is  used  to 
obtain  a  lower  bound  on  reliability  when  the  number  of  non-K-vertices  is  too  large. 
We  demonstrate  that  this  lower  bound  will  usually  be  "good.’’ 

Chapter  5  describes  an  algorithm  to  compute  R^c(D)  or  find  a  lower  bound  on 
R*k(D)<  and  its  implementation.  The  program  is  coded  in  FORTRAN  with  one  main 
routine  and  18  subroutines  totaling  about  1600  lines  of  code. 

Chapter  6  gives  computational  experience  for  the  algorithm  run  on  an  IBM 
3033AP  system  /370  with  the  FORTRAN  H  (Extended)  compiler. 


Chapter  7  is  a  conclusion  and  gives  suggestions  for  further  research. 


ii.  DErramoNs  and  network  rkijahtuty-  theory 


In  this  chapter,  we  give  the  definitions  and  notation  used  throughout  this 
thesis,  and  introduce  basic  network  reliability  theory. 

A.  DEFINITIONS  AND  NOTATION 

A  directed  network  D=(  V.E)  comprises  two  finite  sets:  V  is  the  set  of  vertices 
and  E  is  the  set  of  edges.  Each  edge  eeE  corresponds  to  an  ordered  pair  of 
vertices  in  the  directed  network.  That  is,  e=(u,v )  where  u,V€V.  The  vertex  u  is 
called  the  vertex  of  edge  a  and  vertex  v  is  called  the  head  vertex  of  edge  e  . 
For  any  two  vertices  u  and  v  in  D,  if  there  exists  an  edge  a  ={u,v),  then  u  and  v 
are  said  to  be  adjacent  and  edge  a  is  incident  into  v  and  incident  out  of  of  u  The 
indegree  of  a  vertex  v,  denoted  indeg(v).  is  the  number  of  edges  incident  into  v 
and  the  outdegree  of  a  vertex  v,  denoted  outdeg(v),  is  the  number  of  edges 
incident  out  of  v .  Two  vertices  are  connected  (or  "communicate”)  if  there  exists  a 
sequence  of  vertices  and  edges  of  the  form  v0,  (v0vj).  «i.  (t'i,V2),...((vm_l,i;IR)  vm. 
This  sequence  is  a  path  of  length  m.  If  the  v0=vm,  the  path  is  a  cycle. 

Consider  a  undirected  graph  G(V,E)  formed  by  ignoring  the  ordering  of  the 
edges  in  a  directed  graph  D{V,E).  G  is  called  the  underlying  graph  of  D.  D  is 
connected  if  its  underlying  graph  G  is  connected.  G  is  connected  if  there  is  an 
undirected  path  between  all  pairs  of  vertices.  The  vertex  connectivity  of  D  is  the 
minimum  number  of  vertices  which  must  be  deleted  from  G  (along  with  adjacent 
edges,  of  course)  in  order  to  disconnect  G  or  to  create  a  trivial  graph  with  only  one 
vertex  D  is  separable  if  its  vertex  connectivity  is  one.  Otherwise  D  is 
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nrai —parable.  In  a  separable  graph,  a  vertex  whose  deletion  disconnects  the  graph 

is  called  a  cut-vertex 

Merging  a  vertex  v  into  another  vertex  w  means  removing  all  edges  between  v 
and  w  and  forming  a  single  vertex  v'  by  coalescing  v  and  w  such  that  any  edge 
incident  into  (out  of)  v  or  w  is  now  incident  into  (out  of)  v'.  Merging  v  and  w  is 
sometimes  also  referred  to  as  "identifying”  v  with  u». 

The  SKT  reliability  of  a  directed  network  D,  denoted  •/?**( Z?).  is  the  probability 
that  all  K- vertices  in  D  are  connected  by  path  of  working  edges  and  vertices  from 
source  vertex  s.  The  probability  that  edge  et  works  is  =  l-?<.  For  now.  wt 
assume  that  all  vertices  work  with  probability  one. 


B.  NETWORK  THEORY  AND  SURVEY 


We  briefly  discuss  the  complexity  of  the  SKT  reliability  problem  in 
terms  of  the  theory  of.  NP-completeness  and  in  terms  of  practical  computation. 
More  detail  can  be  found  in  Garey  and  Johnson  [1979],  Ball  and  Provan[1981]. 

Algorithms  to  solve  various  problems  can  be  broadly  classified  into  two 
categories,  polynomial-time  and  exponential-time  algorithms.  An  algorithm  is  a 
polynomial-time  algorithm  if  for  a  problem  of  size  n ,  its  running  time  is  bounded 
by  a  polynomial  in  n .  Any  algorithm  that  is  not  a  polynomial-time  algorithm  is  an 
exponential-time  algorithm.  In  combinatorics,  a  decision  problem  P  is  said  to 
belong  to  the  class  NP  if  given  a  tentative  solution,  the  validity  of  the  solution  can 
be  checked  in  polynomial-time.  P  is  "NP-complete"  if  it  is  equivalent  to  the 
satisfiability  problem  of  a  conjunctive  normal  form  boolean  expression.  A  problem 
is  NP-hard  if  it  is  at  least  as  hard  as  an  NP-compiete  problem. 


A  problem  is  #P-complete  (number-P  complete)  if  its  solution  is  equivalent  to 
counting  (not  enumerating)  the  number  of  solutions  to  an  NP-complete  problem. 
Any  IP-complete  problem  is  NP-hard. 

2.  Enahlam  CompLexily 

Bail  and  Provan  [1981]  show  that  the  ST  reliability  problem  in  an  acyclic 
directed  network  is  a  #P-complete  problem.  This  implies  that  the  SKT  problem  is 
also  IP-complete  since  problem  instances  restricted  to  \K\-Z  are  IP-complete. 
Even  finding  an.  £-approximation  to  SKT  reliability  is  an  IP-complete  problem. 

In  practice,  these  complexity  results  imply  that  a  general  problem  of 
large  size  cannot  be  solved  analytically  in  a  reasonable  amount  of  time.  In  fact, 
the  time  taken  grows  exponentially  with  size  of  the  network.  As  will  be  seen, 
however,  practical  problems  of  fairly  large  size  cam  be  hemdled  with  the  techniques 
developed  here.  Furthermore,  even  larger  problems  cam  be  handled  if  one  is 
willing  to  accept  a  "good”  lower  bound  on  exact  reliability. 

3.  Ppriiirtir.n 

Many  reduction  schemes  have  been  developed  to  help  solve  various 
network  reliability  problems.  These  reductions  reduce  the  size  of  the  network  in 
polynomial  time  while  preserving  reliability,  i.e.,  they  create  a  new.  smaller, 
network  5*  from  the  original  network  G  ( G  may  directed  or  undirected  and  any 
number  of  definitions  of  reliability  may  be  used.)  such  that  R{G)-M  R(G')  where  M 
is  a  known  constamt  factor.  In  this  way,  a  polynomial  amount  cf  work  is  expended 
in  order  to  save,  it  is  hoped,  an  exponential  amount  of  work  during  later,  more 
general  computations.  In  some  special,  but  often  practical  cases,  it  is  possible  to 
completely  reduce  a  network  and  effectively  compute  its  reliability  in  polynomial 
time. 

The  reductions  describe  below  give  a  flavor  for  the  different  types  of 
reductions  which  have  been  developed.  Parallel  and  series  reductions  are  well- 
known  schemes  so  they  will  not  be  discussed.  The  reductions  described  assume 
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that  vertices  are  completely  reliable.  However,  most  can  be  extended  to  handle 
unreliable  vertices. 

(1)  Polygon-To-Chain  Reductions:  A  chain  C  in  an  undirected  graph  G  is  defined  as 
a  alternating  sequence  of  vertices  and  edges  beginning  and  ending  with  vertices 
such  that  all  internal  vertices  have  degree  2.  Two  parallel  chains  form  a  polygon 
and  this  polygon  cam  be  replaced  by  a  chain  in  a  reduction  that  preserves 
reliability  in  the  K-terminal  reliability  problem,  the  analog  to  the  SKT  reliability 
problem  in  an  undirected  network.  This  reduction  and  its  implications  are 
developed  by  Wood  [1982]. 

(2)  2-Neighbor  Vertex  Reduction:  For  the  SKT  reliability  problem,  any  vertex  u 
that  has  exactly  two  neighbors  v  and  w  can  be  eliminated  from  the  network  if 
uf/  K  or  \u,v,w  ]eK  and  u*s  . 

(3)  Parallel-Chain  Reduction:  Suppose  2Ty  is  the  set  of  all  edges  in  a  directed 
network  D  whose  endpoints  are  and  In  D,  a  chain  C.  is  an  alternating 
sequence  of  vertices  and  non-empty  sets  of  edges  where  the  internal  vertices  are 
2-neighbor  vertices  in  D.  Two  chains  C,  and  Cz  with  identical  endpoints  are 
parallel  chains.  Parallel  chains  can  be  reduced  to  a  single  chain  in  the  SKT 
reliability  problem.  Reductions  (2)  and  (3)  are  developed  by  Agrawal  [19B2]. 

4.  nenompositinn 

Sometimes  it  is  possible  to  divide  a  given  graph  into  a  number  of 
subgraphs,  such  that  each  subgraph  can  be  analyzed  separately  and  the  results 
combined  to  obtain  the  reliability  of  the  whole  graph.  We  briefly  describe  some  of 
these  methods  based  on  vertex  connectivity. 

(1)  2-Connectad  Graph  Decomposition;  Let  G  be  a  2-connected  graph  with  a  pair  of 
vertices  \u,v\  whose  deletion  disconnects  the  graph.  These  vertices  are  called  a 
"separating  pair.”  This  separating  pair  partitions  G  into  two  subgraphs  G\  and  G z 
such  that  =G.  G\ C\Gz  =\u.v\,  E and  It  has  been  shown  by  Wood 

[1982]  that,  in  the  K-terminal  reliability  problem,  each  subgraph  defines  a  chain 


through  a  set  of  reliability  computations  on  the  subgraphs.  The  subgraphs  may  be 
replaced  by  the  chains  so  as  to  preserve  reliability.  Agrawal  [1982]  has  proven  the 
analogous  decomposition  scheme  for  the  SKT  reliability  problem  on  directed 
networks.  (2)  Strong-connected  Decomposition;  A  strong  component  of  a  directed 
graph  D  is  a  maximal  subgraph  of  D  such  that  there  exists  a  directed  path  from 
every  vertex  to  every  other  vertex  in  the  subgraph  When  computing  R3x{D),  each 
strong  component  can  be  considered  separately.  Computing  each  component 
R,jc{Di)  and  multipling  /?**•(  A)  is  whole  graph  reliability  /?,*(£).  This 
decomposition  has  been  developed  by  Shogan  [1973]  and  Agrawal  [1962]. 

Other  decomposition  schemes  exist.  Decomposition  in  3-connected 
graphs  is  often  possible.  Rosenthal  [1974].  Decomposition  in  1-connected  graphs 
will  be  discussed  in  chapter  3. 


Any  algorithm  to  solve  a  general  network  reliability  problem  has  a  running 
time  which  is  exponential  in  the  size  of  problem.  If  by  some  method,  the  size  of  the 
problem  can  be  reduced,  immense  computational  savings  can  be  made.  Therefore, 
in  this  chapter,  we  introduce  reduction  and  decomposition  schemes  which  will 
often  reduce  the  size  of  the  problem  while  always  preserving  reliability.  Some  of 
these  techniques  are  specific  to  the  SKT  reliability  problem  in  acyclic  directed 
networks  while  others  are  of  more  general  applicability. 

Reduction  schemes  reduce  a  parameter  of  the  size  of  the  problem,  usually 
edges  or  vertices,  in  order  to  decrease  computational  complexity  while  preserving 
network  reliability.  On  the  other  hand,  decomposition  schemes  split  the  whole 
graph  into  subgraphs,  computes  subgraphs  reliabilities  separately  and  then 
combines  these  reliabilities  in  some  way  to  compute  overall  network  reliability. 

A.  REDUCTIONS 

In  this  section,  we  discuss  various  reduction  schemes  which  reduce  the  size  of 
D  and.  consequently,  reduce  the  complexity  of  computing  RS^{D).  Under  any 
reduction,  reliability  of  the  network  remains  invariant  up  to  a  known,  constant 
factor,  i.e.,  =  M'Rtfc(D')  where  M  is  the  known  factor,  fC  is  the  new  set  of 

K-vertices,  and  D'  is  the  reduced  network.  In  the  following  illustrations,  all  vertices 
are  represented  as  circles  with  K-vertices  being  darkened  or  shaded. 

The  parallel  and  series  reductions  are  standard,  i.e.,  well-known  and  widely 
used.  The  three  reduction  schemes, 


‘'indegree-one  K-vertex  contraction",  "non-K-vertex  deletion"  and  "neck  vertex  to 
K-vertex  reduction"  are  newly  developed  reductions. 

1.  Parallel  Edges  Rednrtinn 

In  a  directed  network,  two  parallel  edges  ea={u,v)  and  e6=(u,v)  may 
be  replaced  a  single  edge  ae  -(u,v)  while  preserving  reliability  if  pc  =  pa+pb  ~PaPb 
and  M~l. 

2.  Series  Edges  Rediirfinn 

In  a  directed  network,  a  pair  of  edges  e„=(u,tO  and  e6=(v,uj)  such  that 
indeg(i»)=outdeg(u)=l  are  series  edges.  If  K,  edges  e„  and  eb  can  be  replaced 
by  single  edge  ec=(u,ui)  while  preserving  reliability  if  pe  =  papb  and  M=l. 

3.  Indegree-One  K-veri-ev  Printr  art-inn 

Let  ei=(u,v)  be  an  edge  in  D  such  that  veK  and  indeg(v)=l.  In  the  SKT 
reliability  problem,  communication  from  s  to  vertex  veK  can  only  only  occur  via 
edge  e<  so  this  edge  must  work  if  the  network  is  to  work.  Consequently,  edge 
can  be  contracted  such  that  R^lD)  -  M'RU(D')  and  M=pt 


where  D’  =  j  V  -u  -v  +  w,E  -  e*  j 
w  =  u{jv 


1C  = 


K-v+w,  itutf  K 
K-v-u+vj,  if  ueK. 


The  validity  of  this  reduction  is  easily  shown.  Let  ei(u,v)  be  the  edge 
described  above  in  the  acyclic  network  D.  Let  denote  the  event  that  e*  is 
working,  and  let  Fi  denote  the  complementary  event.  Since  Rsx{D)  is  just  a 
probability,  the  rules  of  conditional  probability  can  be  applied  to  obtain 


R*(D)  =jnR*(D\Fi)+<hB4c(D\Fi) 

=  PiR*K{D\Ji)  since  R9/dD  j  Ft)=0 


=  PiR*Jc(£') 

since  no  new  paths  are  created  by  contraction. 

4.  Nfm-y-Vartpy  Di»lfttinn 

In  the  SKT  reliability  problem,  an  outdegree  zero  non-K-vertex  can  be 
deleted  from  the  network  without  changing  network  reliability  since  such  a  vertex 
is  irrelevant.  In  this  case  R,x(D)  =  M  RtK(D-v)  and  M-l. 

5.  Mark-  VgrtftV  to  K-rartpy  RpHnrHnn 

Suppose  u  and  v  are  two  vertices  in  D  such  that  ueK,  K  and  all 
paths  from  s  to  u  include  v.  v  is  called  a  neck  vertex.  It  may  or  may  not  be  a 
cut-vertex.  Since  s  must  communicate  with  v  before  it  communicates  with  u , 
R,x(D)  =  U  Rn(K+v)(D)  where  M-l.  Thus,  the  neck  vertex  to  K-vertex  reduction 
simply  changes  K  to  K+v.  See  Figures  3.1a  and  3.1b.  (All  figures  are  given  in 
Appendix  A.) 

The  neck  vertex  to  K-vertex  is  termed  a  "reduction"  because  it  reduces 
the  parameter  |  V-K |  in  D.  As  will  be  seen  in  chapter  4.  the  general  algorithm  for 
computing  JR,k(D)  runs  in  time  proportional  to  Z^V~K]  so  this  is  a  reduction 
according  to  our  definition.  Furthermore,  the  neck  vertex  to  K-vertex  reduction 
will  often  allow  other  reductions  or  decomposition  to  take  place. 

6.  Other  Reriuntinas 

Numerous  other  reductions  exist  as  exemplified  in  chapter  2.  section 
B3.  However,  none  of  these  reductions  has  been  implemented  because  of  their 
complexity,  and  therefore,  no  other  reductions  will  be  described  here. 

B.  DECOMPOSITION 

Sometimes  it  is  possible  to  decompose  a  given  network  D  into  a  number  of 
subnetworks.  Each  subnetwork  can  then  be  analyzed  separately  and  the  results 
combined  in  some  way  to  obtain  overall  network  reliability. 
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Decomposition  can  result  in  substantial  savings  because,  in  general,  the  time 
to  analyze  the  whole  system  is  exponential  in  size  of  the  system.  Effort  needed  to 
compute  Rtf(D)  using  decomposition  is  the  sum  of  efforts  needed  for  each 
subproblem  plus  the  effort,  normally  polynomial  in  the  size  of  D ,  in  performing  the 
decomposition  and  recombining  the  result.  Without  decomposition,  the  necessary 
work  is  roughLy  the  product  of  the  effort  needed  to  solve  the  individual  problems. 


A  connected,  undirected,  network  G{  V,E)  Hm  a  cut-vertex  v  if  there 
exist  components  (subgraphs)  Ct  and  Ca  such  that  Cx  n  C2=v  and  C\ijC2=  G  .  Any 
communication  from  one  component  to  the  other  must  pass  through  v .  In  this 
case  we  also  say  that  v  separates  Cx  and  C2 •  A  network  which  has  a  cut-vertex  is 
"separable"  and  one  which  has  none  is  called  "nonseparable."  A  directed  network  D 
is  separable  if  its  underlying  graph  G  is  separable. 

Let  D  have  nonseparable  components  CX.C2.  .  .  .  C^.  Any  component  Q 
which  has  only  one  cut-vertex  is  called  a  pendant  component.  If  a  pendant 
component,  Q.  contains  no  K-vertices  except  possibly  its  cut-vertex,  then  it  is 
completely  irrelevant  with  respect  to  reliability  since  none  of  its  edges  are  in  any 
path  from  the  source  to  a  K-vertex.  Thus,  this  component,  excluding  its  cut-vertex, 
can  be  deleted  from  the  network  while  preserving  network  reliability  in  the  SKT 
reliability  problem.  This  type  of  deletion  may  be  recursively  applied  ur til  only 
relevant  components  remain. 

Cut -vertices  are  not  efficiently  found  in  a  directed  network  if  each  edge 
is  represented  in  one  direction,  e.g.,  if  all  edges  incident  into  each  vertex  are 
stored  contiguously.  However,  if  edges  are  stored  both  by  tail  vertex  and  head 
vertex,  then  this  is  exactly  the  representation  of  the  underlying  undirected 
network.  Thus,  an  efficient  algorithm  fcr  finding  cut-vertices  in  an  undirected 
graph  may  be  used  to  find  cut-vertices  in  D.  A  slight  modification  to  the 


DFS  (depth-first  search)  algorithm  due  to  Tarjan  [1972]  cam  find  all  cut-vertices  in 
0(  | iS*  | )  time. 

Once  all  cut-vertices  have  been  found  and  all  irrelevant  components 
removed,  any  cut-vertices  remaining  which  are  not  K-vertices  are  neck  vertices 
and  may  be  added  to  the  set  K.  It  is  fairly  obvious  that  each  component  must  work 
by  itself  if  D  is  to  work.  Since  edge  reliabilities  are  independent  it  follows  that 
Rag(D)-Rtx(Di)  R^c{Dz)  •••  An)  where  the  A  are  the  nonseparable, 
relevant  components  of  D.  This  is  identically  equal  to  a  network  which  has  the 
same  components  as  D  but  whose  only  cut-vertex  is  the  source.  A  more  rigorous 
demonstration  of  this  identity  can  be  made  by  applying  the  decomposition  scheme 
described  next,  "moving  edges  to  source,"  to  all  edges  incident  out  of  cut-vertices. 
Figure  3.2  demonstrates  the  results  of  cut-vertex  decomposition  applied  to  a 
network  D  with  five  components,  two  of  which  are  irrelevant.  In  this  case, 
R.jc(D)  =  *«*(  A  ,)•/?.*■  (As)  ( As). 

2.  Moving  Friges  In  finnrnf) 

In  this  section,  we  discuss  a  technique  which  does  not  fit  neatly  under 
the  rubric  of  decomposition  or  reduction  since  it  only  moves  edges  from  one  place 
to  another  in  the  network.  However,  by  moving  certain  edge  tails  to  the  source  s , 
new  parallel  reductions,  series  reductions  or  contractions  may  be  made  possible 
along  with  cut-vertex  decompositions.  Consider  an  edge  e=(u,i/)  such  that  ueK. 
Edge  e  can  be  made  incident  out  of  any  other  K-vertex  u’zK  without  affecting 
reliability  as  long  as  no  cycles  are  created.  A  simple  conditioning  argument  can  be 
used  to  show  this,  Agrawai  [1982].  The  source  is  the  obvious  candidate  to  receive 
the  edge  tails  since  it  has  no  edges  incident  into  it  and  therefore  no  cycles  can 
possibly  be  created.  Figure  3.3  shows  an  example  where  this  operation  creates  a 
cut-vertex  and  allows  a  parallel  reduction. 


IV. 


.►»ri  wl*i  » i  c*lHir+>.  **  *  w > jciw  iuskm 


This  chapter  gives  a  general  method  for  computing  SKT  reliability  in  an  acyclic 
directed  network.  The  method  is  based  on  repeated  computation  of  SAT  reliability 
of  certain  subgraphs  of  D  and,  consequently,  we  first  discuss  SAT  reliability.  After 
deriving  a  method  for  computing  SKT  reliability,  a  simple  lower  bound  becomes 
apparent  and  several  heuristic  arguments  are  given  indicating  that  the  bound's 
accuracy  should  be  good.  Actual  computational  examples  are  reserved  for  the 
next  chapter. 

A.  SOURCE-TO-ALL  TERMINAL  RELIABILITY 

Let  D(V,E )  be  an  acyclic  directed  network  and  let  g<  =  l-Pi  be  the  failure 
probability  for  any  edge  e^E  and  let  Ej  be  the  set  of  edges  incident  into  Vj.  Also, 
let  the  vertices  be  numbered  such  that  •  •  •  ,vn  j  is  an  an  acyclic  ordering  of 

V.  Then.  SAT  reliability  may  be  expressed  as 

R,v(D)  -  P{s  communicates  with  ail  ve.V-vn  (4.1) 

and  at  least  one  edge  into  vn  works  ) 
since  there  no  edges  incident  out  of  vn  by  the  acyclic  ordering 
=  /?,(v-yn)(Z?-a<n)  1-  by  independence 

®  n  1_  FI 
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Another  way  to  see  this  is  to  use  the  fundamental  topological  property  that 
Rtv{D)  is  unchanged  if  any  edge  (u,v)  with  ueK  is  replaced  by  an  edge  ( w,v )  with 
the  same  reliability  and  we/C  has  a  lower  number  than  u  in  an  acyclic  ordering  of 
D.  Consequently,  in  the  SAT  problem,  we  may  move  all  edges  so  that  they  come 
out  of  the  source  just  like  the  "moving  edge  to  source"  operation  in  section  3.2. 
Equation  (4.1)  is  then  trivially  true  for  a  network  where  all  edges  come  out  of  the 
source.  (See  Figure  4.1).  The  complexity  of  computing  Ray(D)  is  easily  seen  to  be 
0(  |  E  |)  if  all  edges  are  with  common  head  vertices  are  stored  contiguously. 

B.  SOURCE-TO-K-TERMINAL  RELIABILITY 

Now  we  consider  the  SKT  reliability  problem  in  D.  Let  D(V,E)  be  a  acyclic 
directed  network  with  KqV.  By  simple  state-space  partitioning, 

/?rt(Z?)  ~  P{s  communicates  with  all  ve V) 

+  ^  P{  s  communicates  with  all  ueK— vx  but  not  within) 

+■  £  P{  s  communicates  with  all  ve7— Vi~Vj  but  not  with  or  Vj) 

vlMj*  K 

+  .  .  . 

+  P{  s  communicates  with  all  v  e K  but  not  with  any  v  e  V-K) 

But. 

P{  s  communicates  with  all  veV)  =  Ray(D) 

P{  s  communicates  with  all  veV— vx  but  not  with  vx  )  =  R^y-vjiD).  P{  No  edge  into 
Vi  works ) 

P{  s  communicates  with  all  veV^vi-Vj)=R,^v-vt-vj)(D)-  P{  No  edge  into  vx  or  vj 
works  except  an  edge  e-{vx,v j)  or  e={vj,vx)) 


R*{D)  =  R,v(D)  +  j  [,^i)R^D  “**)] 


+  .  .  . 


where  £*  is  the  set  of  edges  going  into  vit  £y  is  the  set  of  all  edges  going  into  or 
vj  except  edges  going  between  i/*  and  vj.  E^  is  the  set  of  all  edges  going  into  vit  Vj 
orvh  except  edges  going  between  pairs  of  vertices  in  \vi%Vj,vk },  etc. 

The  complexity  of  computing  Rax{D)  is,  via  the  above  formula,  exponential  in 
the  number  of  non-K-vertices.  The  total  number  of  additive  terms  in  equation  (4.2) 

is  2  Q1  V~K^  -  2,K-Jr|  where  C?  is  the  number  of  ways  to  choose  i  elements  out  of  n 

t«0 

elements.  Computing  SKT  reliability  thus  requires  0(|iT|2,K~ir|)  time  using 
equation  (4.2).  Of  course,  if  the  number  of  non-K-vertices  is  fixed  as  |  £*  |  changes, 
computation  is  effectively  0(  |  E  | ). 

C.  A  LOWER  BOUND  ON  SKT  RELIABILITY 

Since  computation  increases  exponentially  in  the  worst  case,  we  may  be 
satisfied  with  an  approximation  to  SKT  reliability  which  requires  less  time  to 
compute  than  exact  reliability.  Since  all  terms  in  equation  (4.2)  are  positive,  any 
subset  of  these  terms  provides  an  approximation  to  SKT  reliability  which  is  a  lower 
bound.  Also,  the  product  formed  by  lower  bounds  on  independent  components  of  a 
network  yields  a  lower  bound  on  overall  network  reliability.  The  subset  of  terms  we 
choose  to  form  -his  lower,  bound  consists  of  terms  enumerated  in  the  order  given  in 
equation  (4.2)  until  some  specified  computational  limit  is  reached.  The  motivation 
for  this  scheme  being  fairly  accurate  is  given  below. 


«r. •' 


Suppose  Pi  =  p  for  all  and  p  Is  close  to  1.  Consider  the  first  few  terms  of 


Equation  (4.2). 

R.v(D)=  njl-g1'®’*1]  ~  1.0 

v#cV*  ' 


,1*1 


Thus,  the  contribution  of  the  terms  to  the  total  reliability  sum  diminishes 
rapidly  unless  there  are  a  significant  number  of  edges  connecting  non-K-vertices. 


Another  reason  that  the  lower  bound  should  often  be  good  is  that  terms  with 
many  vertices  deleted  may  be  identically  zero.  For  example,  consider  an  ST 
reliability  problem  where  the  shortest  path  from  s  tot  has  1  edges  in  it.  This  path 
must  contain  1-1  non-K-vcrtices  and  so  R,(v-vi-vj-...-vk){D— Vi-Vj-...-vk)  s  o  when 

|  V-fa.Vj . vk  s,t  j  |  <1  -1.  Of  course,  other  terms  may  be  identically  zero,  too. 

For  example,  in  any  K-terminal  problem,  the  term  corresponding  to  D  — vt  — will 
be  identically  zero  if  all  paths  from  s  to  some  K-vertex  v  must  include  either  i\  or 

vi- 

Computational  experience  in  Chapter  V  lends  credence  to  the  accuracy  of  this 


lower  bound. 


The  algorithm  presented  here  for  reliability  evaluation  of  directed  acyclic 
networks  is  based  on  the  theory  developed  in  chapter  3  and  chapter  4.  The 
program  is  intended  for  general  use  so  it  was  coded  in  the  widely  available 
language,  FORTRAN.  This  Language  is  available  on  most  computer  systems  and  is 
still  one  of  the  most  popular  languages  among  operations  research  analysts. 

A.  ALGORITHM 

The  objective  in  this  section  is  to  develop  an  efficient  algorithm  for  computing 
SKT  reliability  in  a  directed  acyclic  network.  Initially,  polynomial-time  reductions 
and  decomposition  are  used  to  reduce  network  parameters  and  decompose  the  the 
network  into  subnetworks,  if  possible.  If  reliability  is  not  completely  computed  via 
reliability-preserving  reductions,  the  general  SKT  reliability  expression  is 
computed  as  described  in  chapter  4.  If  a  specified  CPU  time  limit  is  reached  for 
any  component,  the  computation  is  truncated  and  a  lower  bound  on  reliability  is 
produced. 

There  is  a  one  main  algorithm  and  numerous  subroutine  algorithms  in  the 
program.  VTe  discuss  the  main  program  and  several  important  subroutine 
algorithms  for  understanding.  The  structure  of  the  main  program  is  based  on  a 
step  by  step  approach  which  avoids  redundant  work.  Comments  will  be  given  in 
curly  brackets. 


Input:  A  directed  acyclic  network  D  with  vertex  set  V,  |  7!ss2.  edge  set  E ,  \E  '&2, 
and  set  A'CV,  |  K\^Z.  Edge  reliability  p*  for  each  edge  e*  in  E. 

Output:  R,k(D),  exact  reliability,  if  CPU  time  less  than  specified  amount,  otherwise 

a  lower  bound. 

(1)  {Initialize}  M  -  1.0,  bound=. false. 

(2)  Read  all  edge  data  and  create  network  data  structure 

(3)  {Change  neck  vertices  to  K-verticesj  K  -  Neck  ( D,K) 

(4)  {Perform  reductions  for  whole  network}  Reduct  ( D.K.M ) 

(5)  If  remaining  vertices  |  V\  s>  1.  go  to  end. 

(8)  \Di,Dg . A}=  Decamp (D) 

(7)  For  each  subnetwork  A 

(a)  {Initialize}  lowbnd=. false.,  set  time  =  0. 

(b)  {Move  edge  tails  to  source}  For  each  edge  e=(vi,T;;)  with  v^K 
{move  vt  to  s  j  Move(a) 

(c)  Construct  the  data  structure  for  subnetwork  A 

(d)  {Reduce  subnetwork}  Reduct  (A.A.^f) 

(e)  {  Initialize  }  sumpro=0 

(f)  {  Find  SAT  reliability  }  sumpro  =  Rsa(A). 

(g)  Count  the  remaining  non-K-vertices  for  A.  i-e.,  if  vertex  number 
>  0,  counts count-i-1 

(h)  If  counts  1,  \  find  the  terms  in  . I 

1)  {Find  the  next  combination  of  non-K-vertices} 

h'jj.Uj, . VU]  =  CombO' ,*) 

2)  {  Find  for  A,;a  ;*  I  M'  ~  UaSO{vjt.vJg . «*) 

3)  sumpro=sumpro  +  M'-  Rsa(i7-i/y1— v3a- . — i/^) 

4)  If  elapsed  CPU  time  s  time  limit  then  lowbnd=.true.  and 


bound=.true. 


5)  If  lowbnd=. false.  and  all  combination  have  not  been 
enumerated  then  go  to  (8.h.l) 

(i)  M-M  x  sumpro 

(8)  If  bound=. false,  print  "Exact  reliability  is"  M,  otherwise  print  "Lower 
bound  on  reliability  is"  M 

Bid  Haim 

function  Neck  (D,K) 

Input:  D.K 

Output:  New  set  of  IC-vertices 

(  This  routine  finds  all  neck  vertices  and  changes  them  to  K-vertices  ) 

(1)  Put  all  Vi& K  into  Qn 

(2)  k-\K\ 

(3)  While 

(a)  Remove  v  from  Qn 

(b)  Search  all  K-vertices  in  D-Vi  using  Depth-First  Search 

(c)  k'~  number  of  K-vertices  reached 
(c)  If  Jb'-cfc  then  K-K\j\v] 

(3)  F.eturn  (K)- 

Bid  of  Neck 

subroutine  Reduct(27  ,/T.Af) 

Input:  D.K.M 

Output:  Reduced  network  D  with  modified  K  and  M 

(  This  routine  performs  all  non-K-vertex  reductions,  series-parallel 
reductions  and  indegree-one  K-vertex  contractions  ) 

(1)  Construct  Qd  =  \  veV-K  \  outdeg(i/)=0  j 
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(2)  Whiie  ft  *<f> 

(a)  For  each  vertex  v}  Incident  to  vt,  if  outdeg(v;)  =  l  and  v^V-K 
then  add  to  ft 

(b)  Nonkre  ( i», ) 

(3)  Construct  ft  =  {  v^aV-K  |  indeg(vt)=  outdeg(vi)=l  j  and  ft  =  {  v^K  i 

indeg(ut)  =  l  i 

(4)  For  each 

(a)  Check  all  outgoing  edges,  if  e  ,Vj),e  (vitVf)  are  found  then 

Parall(  e) 

(b)  If  parallel  reductions  creates  new  series  or  contraction  vertices, 
put  these  into  ft  or  ft ,  respectively. 

(5)  While  ft*#  and  ft*# 

(a)  If  ft  *#  Seiies(  ft(i)) 

1)  If  series  reductions  create  new  parallel  edges,  do  parallel 
reduction.  Parali(  a  ) 

2)  If  parallel  reductions  create  new  series  or  contraction 
vertices,  put  these  into  ft  or  ft ,  respectively. 

(b)  If  ft  Contra(ft(i)) 

1)  If  the  contraction  creates  a  new  contraction  vertex,  put 
these  into  ft 

2)  If  the  contraction  creates  new  parallel  edges,  do  parallel 
reduction  Parall(3) 

3)  If  parallel  reductions  create  new  series  or  contraction 
vertices,  put  these  into  ft  or  ft ,  respectively 

(8)  Return  ( D,K,M ) 


Bad  Redact 


subroutine  Decomp(A) 

Input:  D 

Output:  Subnetworks  Dl.Dz . Dm 

(1)  While  not  ail  cut-vertices  have  been  found 

(a)  Using  DFS,  find  pendant  component  A  with  cut-vertex  ut 

(b)  If  A  has  no  K-vertices  except  possibly  vx  then  delete  A 

(2)  Return  (Dlt  Dz . An) 

End  De  comp 


After  reading  in  the  data  and  creating  the  network  data  structure,  the 
algorithm  begins  by  changing  any  neck  vertices  to  K-vertices  by  a  cedi  to  Neck  at 
step  (3).  Next,  a  cedi  to  Reduct  is  made  at  step  (4).  Reduct  first  removes  edl 
outdegree  zero  non-K-vertices.  Through  the  queue  mechanism,  it  is  ensured  that 
all  initial  outdegree  zero  non-K-vertices  are  deleted  along  with  any  which  are 
created  as  the  reduction  proceeds.  A  similar  method  is  used  for  the  other 
reductions  schemes. 

The  other  reductions,  indegree-one-K-vertex  contractions  and  parallel  and 
series  reductions,  may  recursively  create  new  reducible  vertices  or  edges.  The 
newly  created  reducible  vertices  are  put  into  a  queue  instead  of  immediately  being 
reduced  to  avoid  redundant  work.  Parallel  edges  are  reduced  immediately, 
however.  The  network  of  Figure  5. 1  is  completely  reduced  by  reductions  alone  and 
illustrates  how  a  reduction  can  create  other  reducible  edges  or  vertices.  A  call  to 
Reduct  would  result  in  the  set  of  reductions  given  below.  This  network  would  be 
completely  reduced  and  its  reliability  completely  calculated  by  these  reductions. 


(1)  Reduce  parallel  edges 

83,84  to  83 

(2)  contract  vertex 

v0 

(3)  series  vertex 

(4)  contract  vertex 

v2 

(5)  reduce  parallel  edges 

8  g,8 g  tO  83 

(6)  series  vertex 

^7 

(8)  series  vertex 

v7 

(7)  reduce  parallel  edges 

e  10, en  to  e  10 

(8)  series  vertex 

(9)  reduce  parallel  edges 

8  io,  e  12  to  e  jo 

(10)  contract  vertex 

^9 

(11)  reduce  parallel  edges 

82,83  to  82 

(12)  contract  vertex 

*3 

(13)  reduce  parallel  edges 

83,87  to  83 

(14)  contract  vertex  vs 

After  reductions  are  performed,  a  network  will  decompose  into 
subcomponents  if  it  contains  any  cut-vertices.  This  is  done  at  step  (6)  of  Main  by  a 
call  to  Decomp.  The  rest  of  the  calculations  are  performed  iteratively  for  each 
component  D%  under  step  (7).  First,  all  edges  with  tail-vertices  in  K  sure  moved  to 
the  source.  This  may  produce  additional  reductions,  so  another  call  to  Redact  is 
made  for  the  subcomponent.  Then,  the  terms  corresponding  to  Equation  4.2  (for 
that  component)  are  produced  until  ail  such  terms  have  been  enumerated  or  the 
CPU  time  limit  is  reached.  Note  that  this  CPU  limit  is  based  on  the  CPU  time  used 
in  calculating  subcomponent  reliability  and  is  not  based  on  total  CPU  time. 

If  the  CPU  time  limit  is  reached  for  any  of  the  subcomponents,  the  final  value 
obtained  by  the  algorithm  is  a  lower  bound  on  reliability.  Otherwise,  the  value 
obtained  is  exact  reliability.  Choosing  a  CPU  time  limit  will  depend  on  computer 
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facilities  available  and  the  network  under  analysis.  Figure  6.5  shows  the  reliability 
increment  as  a  function  of  CPU  time.  The  lower  bound  obtained  using  a  reasonable 
CPU  time  limit  will  not  differ  much  from  the  exact  reliability  for  this  large  scaie 
network  example 

B.  DATA  STRUCTURE 

There  a  many  ways  to  represent  a  network  on  a  digital  computer.  The  most 
straightforward  is  to  use  an  adjacency  matrix:  If  vertices  vt  and  v}  are  adjacent. 
M(i,j)=l.  and  M(i,j)=0  otherwise.  Matrices  of  the  above  type  are,  in  practice, 
extremely  sparse  and  this  representation  is  very  inefficient  in  terms  of  space.  It  is 
often  inefficient  in  terms  of  execution  time,  too.  particularly  when  one  is 
interested  in  operations  that  require  retrieving  information  from  all  vertices 
adjacent  to  a  particular  vertex,  e.g.,  Neck.  Decomp.  Re  duct,  and  Rsa.  In  this 
situation,  since  most  of  the  eiements  of  M(i.j)  are  null,  the  algorithm  will  spend  a 
great  amount  of  time  retrieving  and  comparing  zero  values.  For  a  static  network, 
an  efficient  representation  is  a  packed  matrix.  In  our  algorithm,  the  network  is 
dynamic,  since  it  is  reduced  and  decomposed  many  times  during  execution. 
Therefore,  a  multi-linked  data  structure  (multiple  linked  lists  with  links  to  other 
data  structures  such  as  arrays)  is  used  to  represent  the  network  Every  vertex  has 
a  linked  list  of  adjacent  vertices,  which,  besides  informing  which  vertices  are 
adjacent  to  it.  also  tells  whether  or  not  the  vertex  belongs  to  set  K.  Two  such  lists 
are  kept  for  each  vertex,  an  "adjacent  out"  list  and  an  "adjacent  into”  list.  Of 
course,  each  adjacency  entry  corresponds  to  an  edge.  So,  in  addition,  there  sire 
pointers  indicating  the  addresses  where  information  about  edges  is  kept. 

Figure  5.3  illustrates  this  data  structure.  The  vertex  vx  only  can  send 
communication  to  vertices  v2,  and  v2,  and  vertex  v2  can  receive  the 


communication  from  vertex  u,  and  send  the  communication  to  vertices  v3  and  u4, 
and  so  on  This  structure  is  a  very  efficient  data  structure  for  dynamic  networks. 


C.  ALGORITHMIC  COMPLEXITY 

It  is  important  to  determine  the  complexity  of  the  various  parts  of  the 
algorithm  so  that  overall  complexity  can  be  understood.  In  this  section,  we  analyze 
the  worst-case  complexity  of  the  subroutines,  as  written  and  thus,  find  the 
complexity  of  the  overall  algorithm.  Not  all  of  the  subroutines  have  been  written 
as  efficiently  as  possible  because  of  the  difficulty  in  programming  such  routines 
and  because  of  limited  time.  Furthermore,  for  problems  of  moderate  size,  it  is 
unlikely  that  much  efficiency  would  be  gained  using  more  sophisticated  routines, 
since,  in  practice,  most  of  the  routines  written  seem  to  operate  mere  efficiently 
than  the  worst-case  analysis  indicates. 

There  are  three  different  reduction  schemes  included  in  the  algorithm:  non- 
K-vertex  deletion,  indegree-one  K-vertex  contraction  and  series  and  parallel 
reductions.  However,  the  last  three  of  these  interact  directly. 

The  non-K-vertex  deletion  steps  (1)  and  (2)  in  REDUCT,  requires  OdE"!) 
operations  in  the  worst  case.  For  example,  consider  a  complete  acyclic  graph 
which  has  only  one  K-vertex,  the  source.  Here,  one  vertex  at  a  time  would  be 
deleted  from  the  network,  but,  in  the  process,  every  edge  would  be  examined 
exactly  once  until  a  single  isolated  vertex  remained.  The  initial  queue-building  is  cf 
no  consequence  since  it  is  0(  |  V\). 

Steps  (3)  through  (5)  constitute  the  other  three  reductions.  Queue-building 
again  may  be  disregarded  since  it  is  at  an  0(  |  V| )  operation.  Initial  parallel 
reductions  may  require  0(|£’|)  operations  using  a  bucket  sort  technique  to  check 
for  edges  with  common  end  vertices.  The  central  step  of  the  .algorithm,  step  (5), 
requires  Od/I2)  time.  This  is  true  since  (1)  at  most  |  V\  series  reductions  or 
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contractions  cam  ever  occur  (2)  each  such  reduction  requires  one  0(  |  V\ )  check  for 
a  newly  created  parallel  edge,  and  (3)  once  identified,  actually  carrying  out  the 
parallel  reduction  requires  only  constant  time.  Overall  then,  REDUCT  requires 
0(  |  V\z )  time. 

There  are  two  different  decomposition  schemes  used,  cut-vertex 
decomposition  and  neck-vertex  decomposition.  The  complexity  of  the  cut-vertex 
algorithm  is  OOffl)  since  it  is  based  on  the  0(|£’|)  depth-first  search.  The 
complexity  of  the  neck-vertex  algorithm  will  be  0(  |  V-K\  |£’i)  since  finding  neck- 
vertices  is  based  on  repeated  Q(  \E  |)  searches  on  j D-v  for  each  non-K-vertex  v.  In 
the  worst  case  then,  this  decomposition  requires  0(|  V\  |£|)  time.  However,  if  we 
fix  the  number  of  non-K-vertices  as  described  in  chapter  4,  then  the  complexity  of 
neck-vertex  algorithm  will  effectively  be  0(  |  E  | ). 

It  follows  from  the  above  discussion  that  the  totad  complexity  of  the  reduction 
and  decomposition  algorithms  is  0(lT/|l£’|)  normally,  and  Od^l2)  when  non-K- 
vertices  are  fixed. 

After  reductions,  computing  SKT  reliability  requires  0(  j  E 1 2:  time.  This  is 

obviously  true  since  2 1  combinations  of  non-K-vertices  must  be  produced  and 

an  0(  |  E  |)  SAT  computation  made  for  each  combination.  When  the  number  of  non- 
K-vertices  is  fixed,  the  computation  becomes  0(  |  £*  | ).  The  entire  algorithm 
therefore  requires  0((£’j2,v'flr')  time  in  general  but  only  0(|  V\z)  time  with  non-K- 
vertices  fixed. 
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In  order  to  test  the  efficiency  of  the  SKT  reliability  algorithm,  we  used  two 
kinds  of  networks,  "complete"  acyclic  networks,  and  acyclic  "street"  networks.  A 
complete  acyclic  network  is  a  graph  D(V,E)  such  that  V=\vl,v2.  •  •  •  ,vn]  and 
e=(vitvj)e£  if  and  only  if  i<j.  Any  such  network  with  n  vertices  has  exactly 

M*-1)  edges.  A  network  is  called  anmxn  "street”  network  if  (1)  its  vertices  can 

be  laid  out  in  an  m,  rows  by  n  columns  rectangular  grid/  (2)  for  all  vertices  except 

> 

those  in  the  last  column  there  is  an  edge  incident  out  each  vertex  directed  to  the 
right,  and  (3)  for  all  vertices  except  those  in  the  first  row.  there  is  an  edge  incident 
out  of  each  vertex  directed  upward.  When  n=m  case  the  street  network  is  called 
"square"  (See  Figure  8.1). 

In  the  complete  network  we  assigned  vx  to  be  the  source  and  Gxed  vn  as  a  K- 
vertex.  This  insures  that  all  edges  and  vertices  are  relevant.  Other  X-vertices  were 
chosen  randomly  as  required.  In  the  street  network,  ail  comer  vertices  were  fixed 
as  K-vertices,  the  lower  left-hand  corner  vertex  assigned  as  the  source  and  all 
other  K-vertices  chosen  randomly. 

A.  SOURCE-TO- ALL-TERMINAL  RELIABILITY  COMPUTATION 

As  described  in  chapter  4,  the  complexity  of  computing  SAT  reliability  is 
0(|2f  |).  Since  repeated  computation  of  SAT  reliability  is  the  key  to  computing 
general  SKT  reliability,  it  is  interesting  to  see  how  much  CPU  time  is  required  for  a 
single  SAT  computation.  Figure  8.2  shows  the  CPU  time  used  as  a  function  of  the 
number  of  vertices  in  the  two  types  of  test  networks.  |  V\  is  increased  from  15  to 
80  for  the  complete  network  and  from  9  to  144  for  the  square  street  network.  Of 
course,  as  a  function  of  \V\,  |  £  |  is  quadratic  for  the  complete  network  ranging 


from  105  to  3160.  \E\  increases  linearly  as  a  function  of  |  V\  for  the  street 

network  ranging  from  12  to  284.  From  the  figure,  it  is  clear  that  SAT  reliability  can 
be  computed  very  rapidly. 

B.  SOURCE-TO-K-TERMINAL  COMPUTATION 

Here  we  investigate  the  work  required  by  the  algorithm  for  computing  general 
SKT  reliability  and  the  efficiency  gained  by  using  reductions  and  decomposition. 
Initial  reductions  and  decomposition  require  0(  |  V\  |  E  | )  time.  After  reductions  and 
decomposition,  the  complexity  of  computing  SKT  reliability  is 
0(maxH2|Vl_A*'  I  |£i  i  j)  where  C  is  the  set  of  separable  components  of  D  and  where 

Vi-Ki  is  the  set  of  non-K-vertices  in  component  i.  Consequently,  we  expect  an 
exponential  growth  in  the  worst  case  for  SKT  computations.  Figure  6.3  shows  the 
growth  in  CPU  time  as  the  number  of  non-K-vertices  is  increased,  with  and  without 
reductions  and  decomposition. 

For  the  complete  network,  the  number  of  non-K-vertices  ranged  from  3  to  14 
out  of  a  total  20  vertices.  The  reductions  and  decomposition  do  not  reduce 
computation  much  here  since  the  complete  network  is  so  densely  connected.  On 
the  other  hand,  computation  times  are  significantly  reduced  when  reductions  and 
decomposition  are  used  in  the  street  network,  primarily  because  many  non-K- 
vertices  are  changed  to  K-vertices.  Notice  that  the  computation  time  required  by 
the  street  network  without  reductions  and  decomposition  is  almost  the  same  as  for 
the  complete  network.  This  shows  how  complexity  is  dependent  mostly  on  the 
number  of  non-K-vertices. 

C.  COMPUTATION  WITH  NON-K-VERTICES  FIXED 

If  the  number  of  non-K-vertices  is  fixed,  i.e.  \V-K\=c  for  some  constant  c, 
then  the  complexity  of  computing  SKT  reliability  is  0(|  V!z)+0(c  |iT|)  -  0(2*  ,E  i) 
which  is  0( )  K|  ®).  Thus,  the  complexity  is  polynomial  instead  of  exponential.  Of 


course,  if  reductions  and  decomposition  are  not  used,  the  complexity  becomes 
0(  |  £*  | )  which  is  better  than  0(|F|2).  However,  in  practice,  the  reductions  and 
decomposition  are  very  efficient  and  do  not  seem  to  increase  computation  tune. 
Figure  6.4  shows  the  increase  in  computation  for  both  test  networks  as  |  K]  is 
increased. 

D.  LOWER  BOUND  COMPUTATION 

As  described  in  chapters  4  and  5.  a  lower  bound  on  SKT  reliability  is  obtained 
when  the  normal  SKT  reliability  computation  is  truncated.  Two  heuristic 
arguments  were  given  in  chapter  4  indicating  that  this  lower  bound  should  be  fairly 
good.  Figure  6.5  shows,  for  actual  computations,  how  the  lower  bound  approaches 
the  exact  network  reliability  very  quickly.  The  test  networks  were  assigned 
different  values  of  Pi  to  show  that  the  accuracy  of  the  bound  is  not  much  affected 
by  Pi-  It  appears  that  a  CPU  time  limit  can  be  used  successfully  for  truncating  SKT 
computation  to  obtain  a  good  lower  bound  when  exact  computation  is  not  feasible. 
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This  thesis  has  explored  the  source-to-K-terminal  reliability  analysis  problem 
in  acyclic  directed  networks.  We  have  shown  that  the  polynomial-time  reductions 
and  decompositions  reduce  an  exponentially  complex  problem  to  a  polynomially 
complex  problem  under  certain  conditions.  Even  when  these  conditions  do  not 
hold,  surprisingly  large  networks  can  be  analyzed.  The  theoretical  and 
computational  results  obtained  here  can  be  used  by  researchers  in  analyzing  the 
reliability  of  directed  acyclic  networks  and,  in  some  cases,  for  the  partial 
reliability  analysis  of  general  directed  networks. 

We  assumed  that  all  vertices  work  perfectly  and  that  edge  reliabilities  are 
constant  while  the  network  system  is  alive.  These  assumptions  may  not  be 
acceptable  in  a  real-world  situation.  It  may  be  necessary  to  treat  vertex 
reliabilities  explicitly  and.  in  some  way,  handle  component  reliabilities  changing  as 
a  function  of  time.  This  will  complicate  reliability  analysis  but  does  not  make 
analysis  infeasible. 

The  extension  to  unreliable  vertices  is  not  hard.  All  K-vertices  must  work  if 
the  system  is  to  work  and  thus,  for  unreliable  K-vertices,  we  obtain  network 
reliability  as  usual  and  then  multiply  by  the  probability  that  all  K-vertices  function. 
Unreliable  non-K-vertices  may  be  handled  in  several  ways.  One  of  the  simplest 
ways  is  to  split  each  unreliable  non-K-vertex  into  two  perfectly  reliable  vertices 
with  an  unreliable  edge  between  the  two  vertices  with  a  failure  probability  equal  to 
the  failure  probability  of  the  original  vertex.  All  edges  incident  into  the  original 
vertex  are  made  Incident  into  the  tail  vertex  of  the  new  edge  and  all  edges  incident 
out  of  the  original  vertex  are  made  incident  out  of  the  new  head  vertex.  However, 


it  is  also  possible  to  redefine  the  reduction  and  decomposition  schemes  and 
redefine  the  general  reliability  expression  to  directly  handle  unreliable  non-K- 
vertices.  For  example,  the  series  reduction  can  include  the  vertex  reliability  so 
that  pe  =  PaPbPi.Vi).  Also,  the  term  in  the  general  reliability  expression  associated 
with  "no  edges  in  ^  work"  can  be  replaced  with  a  term  which  expresses  "no 
edges  into  v*  work  or  v._  fails,  and  no  edges  into  vj  work  or  xtj  fails,  and  ...  except 
edges  going  between  pairs  of  vertices  in  {i»t ,Vj ,  ■  ■  ■  ,vk\.  This  is  not  hard  to 
compute. 

Reliability  changing  as  a  function  of  time  can  be  handled  in  one  of  two  ways. 
The  easiest  way  is  to  just  input  component  reliabilities  at  different  times  and 
compute  network  reliability  at  these  times.  Interpolation  can  be  used  to 
approximate  values  between  the  times  for  which  reliability  is  actually  evaluated. 
Another  possibility  is  to  produce  a  symbolic  reliability  expression  that  can  be 
repeatedly  evaluated  by  a  computer  at  different  times,  thus  avoiding  repeated 
reductions,  decompositions,  etc.  This  might  require  large  amounts  of  storage  but, 
if  reliability  can  be  calculated  at  all  in  a  reasonable  amount  of  time,  it  should  be 
feasible.  Once  the  price  is  paid,  computing  reliability  for  any  given  point  in  time 
should  be  very  efficient. 

As  modern  computer,  power  and  other  network  systems  become  larger  and 
more  complicated,  reliability  becomes  a  more  important  part  of  design  criteria. 
Unfortunately,  it  is  obvious  that  as  networks  become  more  complicated  they  also 
become  more  difficult  to  analyze.  Further  research  should  concentrate  on  ways  to 
streamline  computations  and  to  find  accurate  but  easily  computable  bounds.  The 
method  described  in  chapter  4  for  computing  SKT  reliability  could  be  improved  if 
only  non-zero  terms  were  generated;  An  enumeration  scheme  similar  to  that  of 
Satyanaryana  and  Prabhakar  [  1978]  should  help  in  this  respect.  The  lower  bound 
described  in  this  thesis  is  a  good  start  on  finding  accurate  and  easily  computable 
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APPENDIX  A 


v  a  cut-vertex 


Neck  Vertex  to  K-vertex  Reduction 
<Figure  3.1> 


Completely  Reducible  Network  by  Reduction  Scheme 
< Figure  5.1> 


Multilist  Data  Structure 
<Figure  5.2> 


4x4  Square  "Street"  Network 
<Figure  6.1> 
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